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ABSTRACT 

We investigate the morphology and kinematics of the interstellar medium in the 
environs of the open cluster MrkSO, which includes the Wolf-Rayet star WR 157 and 
a number of early B-type stars. The analysis was performed using radio continuum 
images at 408 and 1420 MHz, and Hi 21cm line data taken from the Canadian Galactic 
Plane Survey, molecular observations of the ^^CO (J = 1 ^ 0) line at 115 GHz 
from the Five College Radio Astronomy Observatory and available mid and far IR 
observations obtained with the MSX and IRAS satellites, respectively. 

This study allowed identification of the radio continuum and molecular counter- 
part of the ring nebula SG 13, while no neutral atomic structure was found to be 
associated. The nebula is also detected in images in the mid and far infrared, showing 
the existence of dust well mixed with the ionized gas. We estimate the main physical 
parameters of the material linked to the nebula. 

The interstellar gas distribution in the environs of Mrk 50 is compatible with a 
stellar wind bubble created by the mass loss from WR157. 

The distribution of young stellar object (YSO) candidates in the region shows 
that stellar formation activity may be present in the molecular shell that encircles the 
ring nebula. 
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1 INTRODUCTION 

Massive stars inject large amounts of energy into the inter- 
stellar medium (ISM) through stellar winds, UV radiation, 
and during their supernova (SN) phase at the end of their 
lives. Therefore, the morphology, dynamics and energetics of 
their interstellar environment are strongly modified by the 
presence of these stars. 

Wolf-Rayet (WR) stars are commonly believed to be 
evolved O-type stars which ha ve almost reached the end 
of their nuclear bur ning phase (|Contil 1 19761 . iMaedeij Il983l . 
Ivan der Hucht 200lh . They are characterized by an intense 
mass fiow with m ass loss rates of ~ (1 — 5) x 10~^ Mq 
yr~^ ICappa et al., 2004) and term inal velocities of 1000- 



3000 km s~^ l|van der Hucht 200 ll and references therein) 
which sweeps up the surrounding interstellar matter. The 
gas shed by the star and the swept-up interstellar material 
are piled-up in expanding shells, called interstellar bubbles 
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(IBs). Stellar wind shocks modify the temperature, presure 
and density of the surrounding ISM. The strong UV photon 
flux (with hv ^513.6 eV) of these stars ionizes the bubbles, 
which are detected in the optical and radio ranges as ring 
nebulae. When the ionization front is trapped in the expand- 
ing IBs, these structures have a neutral outer layer which can 
be detected in the Hi 21-cm line emission. The emission dis- 
tribution of this radio line in the environs of these stars has 
shown the presence of cavities and exp anding shells linked 
to the stars and their optical IBs (e.g. IVasguez et al.l |2005| 
and references therein). Molecular gas related to IBs has also 
been found in a number of cases (|Cappa et al.ll200ll ). 

Stellar formation may be favoured in the cold and 
dense outer shells, follow ing t he mechanism pro posed by 
lElmegreen fc Lada (1977)1 and lElmegreen (2000)1 . In their 
model, massive stars excite an Hll region, that expands and 
sweeps up a shell of shocked cold neutral gas. Eventually, 
the dense shell of cold neutral swept-up gas formed around 
these stars fragments and collapses to produce a new gen- 
eration of stars. Up to now, very few studies dealing with 
stellar formation in the molecular shells of IBs have been 
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carried out (see for example lOevI [l996h . and this point is 
still an open question. 

In this paper we report the results of a multifrequency 
study of the gas distribution around WR 157 and its optical 
ring nebula based on both infrared and radio data. Stellar 
formation activity in the region and its relation to the stellar 
wind bubble are also investigated. 



2 WR 157 AND ITS OPTICAL RING NEBULA 



Sh2-157 (|SharplessllT959l ) consists o f two differe nt regions: 
SG13 (= Simeiz274) at b >-0°35' (|Shain 195ll ) which has 
a claw-like appearance, and SG 14, placed at 6 <-0°35', 
which is diffuse and irregular in shape. SG 13 is associated 
with t he open cluster Markar i an 50 ( Mrk 50) llXurner et al.l 
1 19831 . iLwndstron fc Stenholml Il984 ISmith fc WiUid Il994l ) 
located at {l,b) = (111°21',-0°12'), (RA,Dec [J2000]) = 
(23'' 12"', +60° 29')- The brightest member of this cluster is 
the WR star WR 157 {= HD 219460, WN5+B1II, (/,&) = 
(111°19.8',-0°1 3.80, (RA,Dec [J20 00]) = (23''15'" 12.40^ 
+60°27'01.8"'). Ivan der Hucht 20011 '). 

SG 13 and Mrk 50 are located in the Perseus spiral arm. 
Based on CCD UBV{RI)c photometrv. iBaume et al. 20041 
derived a photo metric distance divirks o = 3.46±0.35 kpc for 
the open cluster. Ivan der Hucht 2OO1I gives a spectrophoto- 
metric distance d = 3.4 kpc for WR 157, based on its asso- 
ciation with the cluster. 

Fig. [1] displays the DSSR image of SG 13. The cross 
indicates the position of WR 157. The ring nebula around 
Mrk 50 is easily identified as a claw-like emission region of ~ 
35' X 40' in size, with the open cluster close to the brightest 
section of SG 13. The diffuse emission at galactic latitudes 
b ^ -0°35' corresponds to SG 14. 

From [O III]5007A, [Nil] 6584A a nd [SII] 6717-H673lA 
observations, iLozinskava et al. (1986)1 found that the LSR 
velocity of SG 13 is in the range -53 to -33 km s"\ This 
result is consistent with studies of bright regions in SG 14 
based on t he H166a RRL (- 43 km s"\ Ped lar 1980j, CO ob - 
servations l|Blitz et al.lll982l ). and Ha data (jGeorgelin 1975h . 

The analytical fit to t he circular galactic rotation model 
bv lBrand fc Blitz f 1993)1 predicts that velocities within the 
range -53 to -33 km s~^are located at kinematical distances 
dk — 3.6 - 5.7 kpc, compatible with the spectrophotomet- 
ric distance to WR 157 and M rk 50. A different kinem at- 
ical distance was published b y iGeorgelin et al. f 1973)1 and 
iGeorgelin fc Georgelin f 1976)1 . who found d = 2.5±0.4 kpc. 
In what follows, we adopt a distance d = 3.7±1.2 kpc for 
SG 13. 



3 OBSERVATIONS 

The analysis of the gas distribution was performed using 
observations in the optical, IR and radio ranges. The DSS R 
image of Sh2-157 was obtained from the Skyview web pagc0- 
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Figure 1. DSS R image of SG 13. The cross marks the location 
of WR 157. The grayscale is in arbitrary units. 



3.1 Radio continuum, HI, and CO data 

To analyze the ionized and neutral atomic gas distribution, 
we used Hi observations from the Canadian Galactic Plane 
Survey (CGPS) obtained with the Synthesis Telescope of 
the Dominion Radio Astrophysical Observatory (DRAO) in 
Canada. This telescope performed interferometric observa- 
tions of the 21-cm Hi spectral line, and, simultaneously, con- 
tinuum emission in two bands centered at 1420 MHz and 408 
MHz. The East- West array consists of 7 antennae, 9-m each. 
Single-dish data were routinely incorporated into the inter- 
ferometric images to ensure complete coverage of the emis- 
sion on all angular scales down to the resolution limit. Radio 
continuum observations of the region of SG 13 have synthe- 
sized beams of 3!4 x 3!9 and 58" x 67" at 408 and 1420 MHz, 
respectively. The measured rms image noises are 5 and 1.5 
K at 408 and 1420 MHz, respe ctively. Details about DR AO 
and the CGPS can be found in iLandecker et al. (2000)1 and 
iTavlor et ai]|(2003)l . 

To investigate the neutral hydrogen distribution, we ex- 
tracted a data cube centered at {l,b,v) = (111°8' , -0°23' 
, -40.2 km s"^) from the CGPS. The Hi data have a syn- 
thesized beam of l!l3x0'98, a rms noise of 3 K in bright 
temperature (Ts), and a velocity resolution of 1.3 km s~^ 
with a channel separation of 0.824 km s~^. The Hi images 
were convolved to a 2' x 2' beam size to facilitate the identifi- 
cation of structures. The observed velocities cover the range 
-60 to 163 km s"\ 



The '^CO (J 



0) line data at 115 GHz were 



obtained using the radiotelescope of the Five College Ra- 
dio Astronomy Observatory (FCRAO), in USA. The angu- 
lar resolution is approximat ely 46". Details a bout the CO 
survey are summarized bv .Ridge et al. (2006)1 . 



3.2 IRAS and MSX data 



^ http://skyview.gsfc.nasa.gov/ 



Infrared images at different wavelengths were used to ana- 
lyze the dust distribution in the region. High resolution IR 



The environs of SG 13 3 



"11 




Figure 2. Top panel: Radio continuum image at 408 MHz. The 
grayscale in brighness temperature Tj, corresponds to 110 to 150 
K. The contour lines are 115, 120, 125 130, 140, and 150 K. Middle 
panel: Radio continuum image at 1420 MHz. The grayscale in 
corresponds to 7 to 13 K. The contour lines are 8, 9, 10, 11, 12, 
and 13 K. The cross marks the position of WR 157. Bottom panel: 
Overlay of the optical emission and the 1420 MHz contour lines. 



images (HIRES) at 60 and 100 ^m, were taken from IPAC 
with angular resolutions of ~ 1' and 2'. 

The Spatial Infrared Telescope (SPIRIT III) on board 
the Midcourse Space Experiment (MSX) satellite surveyed 
the entire Galactic Plane in four mid IR bands centered 
at 8.28, 12.13, 14.65, and 21.3 ^m (bands A, C, D and E, 
respectively). We retrieved images in the four bands with 
18'.'4 in angular resolution from IPAC. We converted the 
flux densities to Janskys by using the conversion factor of 
radian ce to flux density given in the MSX Image Server at 
IPAC (|Egan et al.lll999l ). 

To investigate the presence of pre-main sequence stars 
towards the region under study we extracted IR point 
sources from the MSX, 2MASS, and IRAS point source cat- 
alogues. 



4 GAS AND DUST DISTRIBUTION 
4.1 Radio continuum emission 

The radio continuum images at 408 and 1420 MHz are shown 
in the top and middle panels of Fig. [2] The radio emission 
distribution at 1420 MHz correlates strikingly well with the 
optical emission (see the bottom panel of Fig.[2|. The bright- 
est optical emission region coincides with the strongest radio 
emission region in SG 13 at 1420 MHz. In addition to the ion- 
ized branches present at / ~ 111°5' and 111°25', the image 
at 1420 MHz reveals weak radio continuum emission at {I, b) 
~(111°35',-0°15'), coincident with a region of diffuse optical 
emission. Previous high resolution radio continuum studies 
of Sh2-157 were cent ered on SG 14 and did not include the 
northern SG 13 area (|lsrael (1977)1 ). 

The ring nebula around WR 157 can also be identified 
in the radio continuum surveys at 4850 M Hz (Condon et 
al. 1991) and 2700 MHz (|FMrst et al. 1990 I ), taken with the 
Effelsberg 100-m telescope. These images are not included 
in this paper. 

Clearly, the image at 408 MHz is contaminated with 
reduction artifacts from the strong source CasA, lo- 
cated 2° away from Sh2-157. The sources detected at 
(/,6)=(lll°10',-0°30') and («, 6)=(lll°5',-0°22' ), which 
are not evident at 1420 MHz, coincide with the X- 
ray source X IWGA J23138-f6024 and the radio source 
NVSS J231310+601236, respectively. Both are listed in the 
NEClfl database as extragalactic radio sources. These sources 
were not included in the flux density estimate. 

Derived flux densities (S„) at 408, 1420 and 2700 MHz 
are 1.8±0.9, 3.4±0.9 and 3.5±1.0 Jy, respectively. The un- 
certainty in flux density arises in both the rms noise of the 
images and in the estimate of the background emission. 

Uncertainties in the flux densities are too large to allow 
a meaningful determination of the spectral index, however 
the flux densities are consistent with a thermal origin. 

The sources centered at {l,b) = (lll°26',-0°40'), {l,b) 
= (lll°18',--0°40') and {l,b) = (lll°8',-0°46') in the image 
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at 142 MHz were n amed G 111.4-0.7, S 157 A and G 111.2- 
0.8 bv llsrael f 1977)1 . and correspond to SG 14. 

4.2 Infrared emission 

The IRAS (HIRES) image at 60 /im around the open cluster 
is displayed in the top panel of Figure O The middle panel 
displays the MSX band A emission centered at 8.28 fim. The 
7.6 and 8.6 fim features of the policyclic aromatic hydrocar- 
bons (PAHs) strongly contribute to the emission within this 
band. The bottom panel of Fig. [3] shows an overlay of the 
60 fj,Tn (contours) and 8.28 ^m (grayscale) emission distri- 
butions. 

The emission distribution in the far IR presents the 
same morphology as in the optical and radio continuum 
bands. The IR emission at 100 /im, not shown here, displays 
the same distribution as the 60 fim image. The emission in 
the far IR originates in thermal emission from dust. The 
strong correlation between the far IR and the radio contin- 
uum emission suggests that large grains and ionized gas are 
well mixed. 

Strong correlation between the emissions at 60 and 8.28 
fim is clear at {l,b) ~ (111°12', +0°2') and {l,b) ~ (111°5' 
, -0°20'). The 8.28 /xm emission distribution coincides with 
the emission at 60 /xm mainly at the right branch at I ^ 
111°20'. A different behaviour is observed towards (/,&) = 
(111°25',-0°15'), where the bright region at GOfim coincides 
with a low emission region at 8.28 fim. Very probably, the 
dust grains responsible for the emission at 8.28 /im have 
been destroyed in this part of the nebula, which is closest to 
the WR star. 

The strongest IR emission region present at b <-0°30' 
corresponds to SG 14. 

4.3 The emission from the molecular gas 

The CO (1-0) emission distribution of a larger region is 
shown in Fig. |4] The left panel shows the result of integrat- 
ing the emission within the velocity range -58.4 to -43.5 
km s~^ in grayscale and contour lines, while the right panel 
shows an overlay of the contour lines of the left panel and the 
optical emission. A number of CO cloudlets borders the sec- 
tion of the ring nebula toward lower galactic longitudes. The 
brightest optical region, which runs parallel to I — 111° 25', 
appears almost free of molecular material, since only the end 
point near b ~ — 0°5' matches a strong CO cloudlet. The CO 
cloud near {l,b) = (111°25', — 0°32') coincides with the low 
latitude edge of the ring nebula and separates it from SG 14. 
Note that the maximun in the CO emission corresponds to 
an optically faint region. The fragmentary molecular emis- 
sion delineates a roughly circular structure surrounding the 
optical ring nebula, as would be expected where an ionized 
region is surrounded by neutral gas. This is schematically 
indicated by a dashed circle in the left panel of Fig. |3] 

The velocities of the CO gas that surrounds the optical 
nebula span the interval -55 to -43 km s~^. CO associated 
with the region at {l,b) = (111°30',-0°15') was detected be- 
tween -43 and -47 km s^^, while the CO cloud at {l,b) = 
(lll°25',-0°35') is present at ~ -50 km s~^, and the chain 
of GO cloudlets spans the velocity interval -55 to -50 km 
s~^. The region at (i, 6) = (111°35',-|~0°5'), apparently unre- 
lated to SG 13, is detected in the velocity range -58 to -48 




Figure 3. Top panel: IRAS image at 60 fim. The grayscale cor- 
responds to 65-140 MJy/sr. The contour lines are 70, 80, 90, 100, 
110, and 120 MJy/sr. Middle panel: MSX band-A at 8.28 fim. 
The grayscale is 10.5 {7a-) to 40 MJy/sr and the contour lines are 
16, 21, 26, 31 and 36 MJy/sr. Bottom panel: Overlay of the 60 
(in contour lines) and the 8.3/im (grayscale) emissions. 
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km s^^. The positional correlation between the optical and 
CO rim structures and the agreement between the veloc- 
ities of these cloudlets with the corresponding velocity of 
the ionized gas (see Sect. [2]) suggest that the molecular fea- 
ture is the molecular counterpart of SG 13. The clumpy CO 
morphology suggests that most of the CO gas in the region 
has been dissociated by the strong UV photon flux of the 
massive stars in MrkSO. Considering the whole molecular 
velocity range, the expansion velocity of the CO ring is 8±1 
km s~^. 

A large and intense patch of CO emission is present in 
the northern section of the image. With the exception of 
the Hll regioii Sh2-159, located at (l,b) = (lll°36',-t-0°22') 
jlsrael (1977)1 ') . this molecular emission shows little corre- 
lation with the optical emission. The strong CO emission 
projected onto Sh2-159 is very probably associated with the 
Hll region, becau se of the similar CO and Ha velocities 
l|Blitz et al.lll982D . 

A comparison of the IR emission distributions at 8.28 
and 60 nm (Fig. |3} with the molecular one (Fig. |4} reveals a 
clear morphological correlation. The CO clouds encircle the 
IR emission associated with SG 13. 

Adopting a systemic CO velocity of -48 km s~^ for the 
molecular gas related to SG 13, the circular galactic rota- 
tion model by iBrand" fc Blitz f 1993)1 predicts a kinematical 
distance dco ~ 5.0±0.8 kpc. The uncertainty was derived 
adopting a velocity dispersion of 6 km s~^ for the interstellar 
gas. The kinematical distance is compatible with the spec- 
trophotometric distance of Mrk50 and WR157 (see Sect. 
2). 

4.4 The emission of the neutral atomic gas 

Figure [S] exhibits the average Hi emission distribution along 
the line of sight to SG 13. The profile was obtained by aver- 
aging the Hi emission within a box of l°30'xl°45' enclosing 
the Hll region. The more intense peaks are present at and 
-50 km s~^ with brightness temperature Ti, ~ 70 and 90 
K, respectively, enclosing a central region with a brightness 
temperature of about ~ 20 K. A less intense peak is centered 
at ~ -100 km s"^ with Tb ~ 20 K. An analytical fit to the 
circular galactic rotation model by Brand and Blitz (1993) 
predicts that material with velocities of -100 and -50 km 
s~^ should be located at ~12 and ~5 kpc, corresponding 
to the Cygnus and Perseus arms, respectively, while gas at 
about km s~^ corresponds to the Orion arm (local gas). 

Figure |6] displays the superposition of the CO emission 
distribution shown in Fig. [3] (in contour lines) and the Hi 
brightness temperature Tj, distribution within the velocity 
range -48.0 to -33.2 km s~^ in steps of 2.5 km s~^. The 
analysis of the Hi emission distribution within this veloc- 
ity range allows identification of two structures at veloc- 
ities compatible with the radial velocities of the ionized 
and molecular material linked to SG 13. The larger and 
weaker structure, of about ~30' in radius, is centered at 
(Z,6) = (lll°7',-t-0°8') (hereafter named Shell A), while the 
smaller and brighter one, of about 20' in radius, is centered 
at (/, 6) = (110°55',-0°35') (hereafter Shell B). Both shells are 
shown in Fig. |6]as white circles. 

Shell A is clearly identified at ~ -44 km s~^ . The section 
of this Hi structure at 6 ~ -1-0° 20' is associated with intense 
CO emission detected at these galactic latitudes. 



SheU B surrounds part of the Hll region SG 14. The 
section of Shell B at i = 110°47' coincides with relatively 
strong CO emission, while the border at I — 111°20' cor- 
relates with faint CO emission linked to SG 14. However, a 
morphological connection with the IR and radio continuum 
counterparts of SG 13 is not observed. No CO is detected at 
the interface between the two Hi shells. 

The lower right panel of Figure [7] displays the Hi inte- 
grated emission in the velocity interval from -44.0 to -41.5 
km s~^, corresponding to the systemic velocity of Shell A. 
The upper and left panels of Fig.[7]show the {v, I) image cor- 
responding to 6 = 0.0° and the {v, b) image for / = 111°10', 
respectively. These images allow us to estimate the extension 
in velocity of Shell A and are useful to derive its expansion 
velocity. Shell A can be detected within the velocity interval 
from -32 to -50 km s~^, with an expansion velocity of 13±2 
km s-^ Shell B can also be identified on the left panel of 
Fig. [3 although it is less conspicuous that Shell A. Following 
Brand & Blitz (1993) we determine a kinematical distance 
for SheU A of 4.6±0.8 kpc. 

The molecular gas related to the section of SG 13 at 
(/, b) — (111° 10', 0°0') appears projected close to the center 
of Shell A. The distribution of the molecular material re- 
lated to SG 13 correlates neither with Shell A nor with Shell 
B. The analysis of the Hi emission distribution spanning 
the range -48.5 to -51.8 km s~^does not show an obvious 
Hi counterpart to the ionized and molecular gas related to 
SG13. 

In order to obtain a better view of the general HI envi- 
ronment of the two structures, we present in Fig.|8]a slightly 
larger field of view than the one used for the previous Hi im- 
ages. The figure is the average of three CGPS channels be- 
tween -41.86 and -44.34 km s~^, smoothed to a resolution 
of 2'. The figure shows very clearly the shell- like structure 
of feature A. 



5 STAR FORMATION 

In this section, we consider the possibility that star for- 
mation is going on in the expanding shell pushed by the 
stellar winds of the massive stars in Mrk50, particularly 
by WR 157, following the "collect and collapse" model by 
lElmegreen fc Lada (1977)1 . We look for the presence of YSO 
candidates that are projected onto the molecular envelope 
of SG 13 in the MSX, IRAS, and 2MASS point source cat- 
alogues. The searched box area was l°xl° centered on the 

position of the op en cluster. 

We follow the lComeron et ah" (2005)1 criteria in looking 
for YSO candidates in the 2MASS point sources catalogue. 
These criteria discriminate between giant and main sequence 
stars with or without reddening, and sources with IR excess. 
The last ones are the most important sources for our pur- 
pose since their IR fiuxes reveal the presence of circumstellar 
IR emission. According to these criteria, we found around 
20000 2MASS sources, of which only 19 present IR excess. 
Fig. |9] displays the colour-magnitude (CM) diagram of the 
19 sources with IR excess assuming a distance of 3.7±1.2 
kpc. The ZAMS from 03 to B5 type stars is indicated at 
the left of the diagram. Half of the IR excess sources present 
visual absorption larger than 10 mag. Table [l] summarizes 
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Figure 4. Left panel: Integrated CO emission distribution within the velocity range —58.4 to -43.5 km s"'^. The grayscale corresponds 
to 1.5-50 K km s^^. The contour lines are 2.5, 6, 11, 16, 21, and 26 K km s^^. The dashed circle marks the clumpy CO structure. Right 
panel: Overlay of the DSS R image in grayscale and the CO emission distribution in contour lines. 




Figure 5. Average Hi spectrum towards SG 13 
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Figure 6. Superposition of the Hi emission distribution in the line of sight to Mrk50 (grayscale) with the CO integrated emission 
(contours). Each Hi map shows the mean brightness temperature within an interval of ~ 2.5 km s"'^. The grayscale in Tf, of the Hi 
emission is 40 to 100 K for maps (1), (2), (3), and (4); and 70 to 135 K for maps (5) and (6). The cross marks the position of WR 157. 
The white circles delineate Shells A and B detected in Hi. 



the main parameters of these sources, i.e. galactic coordi- 
nates, designation, and J, H, and Ks magnitudes. 

Following the criteria given bv lJunkes et al. (1992)1 . we 
have also found that several IRAS point sources projected 
onto this region are protostellar candidates. The data for 
these sources are compiled in Table [1] which shows the {I, b) 
coordinates, the IRAS name, the fluxes at 12, 25, 60 and 
100 /im of each source, and the FIR luminosity according to 
Chan & Fich (1995). We found that 9 out of the 17 IRAS 
sources are YSO candidates. 

iLumsden et al. (2002)1 derived several criteria to help 
identify massive young stellar objets (MYSOs) in the MSX 
point-source catalogue. The main idea is to discriminate 
among sources with IR excess originating in dust envelopes 
around young and evolved stars, and Hll regions. From mid- 
IR colour-colour diagrams, they found that MYSOs have IR 
fluxes with ratio F21/FS >2 and F14/F12 >1, where Fg, F12, 
Fi4 and F21 are the fluxes at 8.28, 12.13, 14.65 and 21.30 
Hm. For compact Hll regions, these flux ratios are F21/ Fg, > 



2 and F14/F12 < 1. Evolved stars have F21/F8 < 2. On this 
basis, no YSO candidates are found inside the searching box 
area. 

In Fig.[TD]we display the spatial distribution of the YSO 
candidates of Table [1] projected onto the 1420 MHz image 
in grayscale and the CO (1-0) emission distribution in con- 
tours. The crosses correspond to objects with IR excess from 
the 2MASS catalogue and the triangles to IRAS protostellar 
candidates. 

IRAS sources 20, 21, 24, and 26 are projected onto CO 
cloudlets. Sources 25 and 27 appear projected over regions 
lacking molecular emission, while source 23 coincides with 
an ionized filament. On the other hand, 2MASS sources 6, 7, 
11, 18, and 19 are projected over radio continuum filaments, 
while sources 2, 5, 12, 13, and 15 coincide with CO clouds. 
Sources 9, 16, and 17 are projected onto a diffuse ionized 
region. Sources 1, 4, 8, and 22 are related to the Hll region 
SG 14. Finally, sources 14 and 28 are located far away from 
the molecular emission linked to SG 13. 
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Figure 7. Bottom right panel: Hi integrated emission in the velocity interval from -44.0 to -41.5 km s ^, corresponding to the systemic 
velocity of Shells A and B. The grayscale is 80-140 K and the contour lines are 70, 80, 90, 100, 110, 120, 130, and 140 K. Upper panel: 
{v,l) image for b = 0.0° showing brightness temperature T;,. Grayscale: 30 to 90 K. Contours: 30, 35, 40, 45, 50, 60, 70, 80, and 90 K. 
Bottom left panel: {v,b) image for I = 111°9.6' showing the same grayscale and contours as the upper panel. 



We conclude that many YSO candidates are projected 
onto the interstellar bubble envelope. Besides, we believe 
that about half of the 2MASS sources having IR excess 
with j4„ > 10 mag may evolve towards massive stars. Their 
strong visual absorption and their location in the CM di- 
agram favour this suggestion. The existence of a relatively 
large number of 2MASS and IRAS YSO candidates reveals 
the presence of young stellar objets in different evolutionary 
phases. The spatial distribution of the YSO candidates is 
suggestive of the action of the "collect a nd collapse" process 
described bv lElmegreen fc Lada figYTU 



6 DISCUSSION 

6.1 Main physical parameters of the ring nebula 

Table [2] shows the more relevant parameters of the ionized 
and neutral structures linked to MrkSO. 

The electron density rie, the ionized mass Mion, the 
emission measure EM, and the observational excitation pa- 
rameter t/i2d_were_detennm^^ the classical expres- 
sions bv lMezger fc Henderson f 1967)1 from the image at 1420 
MHz. We assumed that SG 13 is a spherical nebula with con- 
stant density. We adopted an electron temperature Te — 10* 
K, and assumed singly ionized He (the derived mass Mion 
was multiplied by 1.27 to take this fact into account). The 
estimated rms electron density and the ionized mass, cor- 
responding to a filling factor / = 1, are listed in Table [S] 
Electron densities for / = 0.1 — 0.2 are also included in 
Table (2] These filling factors were derived by taking into 
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Table 1. YSOs candidates in direction to SG 13 









2 MASS sources 










# 


l[o] 


b[o] 


Designation 


J 

mag 


H 

mag 


mag 






1 


111.41 


-0.56 


23164713-t-6010319 


15.132 


13.265 


11.997 






2 


111.16 


-0.40 


23142490-1-6013421 


11.428 


11.34 


11.209 






3 


111.19 


-0.24 


2314113H-6023109 


8.756 


8.603 


8.399 






4 


111.29 


-0.68 


23161127-1-6001148 


14.002 


12.06 


10.642 






5 


110.91 


-0.22 


23120304+6018385 


13.062 


11.392 


10.147 






6 


111.08 


-0.30 


23133130-f6017273 


13.473 


12.383 


11.363 






7 


111.06 


-0.27 


2313165H-6019056 


11.844 


11.76 


11.63 






8 


111.41 


-0.53 


23164107-1-6011498 


13.713 


12.691 


11.729 






9 


111.61 


-0.33 


23173713-1-6027406 


9.727 


9.727 


9.626 






10 


111.34 


-0.10 


231454234-603425 


12.391 


11.956 


11.582 






11 


111.07 


-0.12 


23125722-f6027223 


11.979 


11.728 


11.505 






12 


111.06 


0.097 


23121250-1-6039495 


10.977 


10.884 


10.739 






13 


111.18 


-0.10 


23134236-1-6031201 


11.932 


11.779 


11.584 






14 


111.72 


0.03 


23172558-f6050436 


7.461 


5.98 


4.826 






15 


111.64 


0.04 


23164590-1-6049490 


15.606 


12.742 


10.971 






16 


111.56 


-0.16 


23164539-t-6035553 


10.929 


10.82 


10.654 






17 


111.37 


0.04 


23144274-t-6043428 


11.194 


10.99 


10.788 






18 


110.99 


-0.51 


23132696-f6003570 


10.932 


10.883 


10.769 






19 


110.90 


-0.41 


23122993-1-6007256 


12.191 


11.742 


11.404 












IRAS 


sources 




















Fluxes 




^FIR 










Ufim [Jy] 


25/im[Jy] 


60/im[Jy] 


100Atm[Jy] 


[Jy] 


20 


110.88 


0.01 


23088+6014 


0.38 


0.41 


6.38 


28.79 


157.92 


21 


111.23 


0.10 


23113+6027 


0.54 


0.83 


8.64 


40.40 


220.55 


22 


111.42 


-0.56 


23146+5954 


7.17 


15. 


275. 


671. 


4324.25 


23 


111.45 


-0.31 


23141+6008 


0.46 


0.21 


8.22 


38.79 


210.13 


24 


111.58 


0.02 


23141+6030 


0.71 


1.05 


9.60 


44.90 


245.51 


25 


111.58 


-0.12 


23146+6022 


0.81 


0.97 


5.88 


32.70 


174.00 


26 


111.63 


0.11 


23143+6036 


0.93 


1.30 


28.39 


116. 


648.31 


27 


111.66 


-0.08 


23151+6026 


0.95 


0.52 


4.40 


23. 


123.88 


28 


111.73 


0.03 


23152+6034 


22.20 


25.10 


147. 


227. 


1800.00 
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GALACTIC LONG. 

Figure 10. Spatial distribution of YSO cantidates in direction to SG 13 over the image at 1420 MHz in grayscale from 7 to 13 K, and 
the ^■^CO image with the following contours: 1.5, 5, 10, 20, 30, and 40 K km s^^. The 2MASS and IRAS sources are indicated by crosses 
and triangles, respectively. 
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Figure 9. CM diagram of 2MASS point sources with IR excess 
in direction to SG 13. 



account that SG 13 is a ring of ~ 12' in radius and 5' in 
thickness, and that only 10%-30% of its surface is covered 
by plasma. 

Based on [Sll] emission lines, iLozinskava et al. figSG)! 
determined electron densities n'^ = 100 — 150 cm~'^ for 
the optical nebula, higher than the values derived from ra- 
dio continuum observations. Following Israel (1978), a dif- 
ferent estimate for the filling factor can be obtained as 
/ = {ne/ne'f'^ = (2-7) x 10"^. Electron densities derived 
from optical line ratios are higher than the ones obtained 
from radio continuum observations. The former corresponds 
to regions with high emissivity and high electron densities, 
while the latter have a large contribution from low electron 
density regions. Very probably, the real / is in between these 
two values. 

f/rad is related to the number of UV photons used to 
ionize the gas and is obtained as C/rad = Rrie^^^, where R 
is the radius of the ionized gas. U, is linked to the number 
of UV photons emitted by the massive stars in the open 
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cluster, which can be obtained from atmosphere models. To 
estimate U, we took into account only WR 157 (WN5 star), 
since t he other star s in Mrk 50 have s pectral types later than 
B3 (|Baume et al. 2004 ). Following S mith et al. (2002)1 . U* ~ 
128 pc cm~^. From Tabled Umd < Ut, indicating that the 
UV photons emitted by the WR star are enough to ionize 
the gas. A large number of UV photons probably escape 
to the ISM through the clumpy molecular envelope around 
SG 13 and/or are absorbed by the associated dust, warming 
and destroying it. 

Bearing in mind that the interstellar dust radiates in the 
far IR, we derived the dust mass Md associated with SG 13 
and the dust colour temperature Td from the emission at 60 
and 100 /im (Draine & Lee 1984). The absorption coefficient 
Xv of the dust can be written as 



Table 2. Physical parameters of SG 13's counterparts 



kg" 



(1) 



where n — 1 
obtained as: 

T4K) 



3 X 1012Hz, 

1.5 is the dust spectral index. Td can be 



95.94 



In [(1.67)3+"^] 



(2) 



5100 and Seo are the IR fluxes at 1 00 and 60 ^m. For 
n = 1 - 1.5, Td ~ 28 to 31 K. Following iHildebrand (1983)1 
the dust mass is 



(3) 



where S^(Td) is the blackbody function. The flux densities 
at 60 and 100 /im, along with the derived dust mass are 
listed in Tabled 

The associated molecular mass of each cloudlet A/h2i 
can be obtained as: 



H2- 



(4) 



We adopted = 2.76 for the molecular weight (assuming 
solar abundances), rriH is the atomic hydrogen mass, Q, is 
the solid angle of the molecular cloudlet and N-n^ is the H2 
column density, which is obtained as 



X X 7co, 



(5) 



where X = 

IIGrenier fc Lebrum 199d . 



1.9x1 0^ 
iDigel et all 



I2OO2I ). 



K kms ^ 

Ico is 

the integrated emission of the CO line. The mean column 
density and the molecular mass are listed in Table [2] 

The gas-to-dust ratio derived taking into account the 
ionized and molecular masses listed in Table [5] is ~ 20, lower 
than the typical value of ~ 100 generally accepted for Hll 
regions. 



6.2 Origin and scenario 

Assuming the CO expansion velocity of 8±1 km listed 
in Table [21 the dynamical age of SG 13, according to wind 
bubble evolutionary models, is td = 0.55J?inoi/i'cxp = 
(1.4 ± 0.4) X 10® yr. Considering the radio continuum, IR 
and molecular counterpart of SG 13, the kinematic energy 
of the interstellar bubble is Sfc=(9.8±6.0)xl0'''^ erg. 

Assuming typical values for the stellar wind of a WN 5 
star, a mass loss rate M = 1.75x10"^ M0 yr~^ and a ter- 
minal velocity Vac = 1520 km s~^ (Smith et al. 2002, Cappa 



Distance adopted 


3.7±1.2 kpc 


Radio continuum 


S270o(Jy) 


3.5±1.0 


Si42o(Jy) 


3.4±0.9 


S408(Jy) 


1.8±0.9 


Spectral index, a 




408 and 1420 MHz 


-l-0.5±0.3 


1420 and 2700 MHz 


-l-0.03±0.1 


Angular radius (') 


12±3 


Linear radius -R(pc) 


15±4 


ne(/=l) 


4±1 


Mio„(/ = 1)(M0) 


3700±600 


ne(/=0.1-0.2)(cm-3) 


15-9 


Mio„(/=O.1-O.2)(M0) 


1100-1500 


t^rad (pc Cm-2) 


45 


U* (pc cm~2) 


~128 


EM(pc cm-'5) 


(2.4±1.7)xl03 


IR 


Sec (Jy) 


~6.9xl03 


Sioo (Jy) 


-2.25x10"* 


Dust color temperature (K) 


28-31 


Dust mass (Mq) 


90±55 


12CO 


(1-0) 


(/, b) centroid of IB 


lll°5',-0°4' 


Velocity range Ai)(km s"^) 


-58.4 to -43.5 


Expansion velocity (km s"-'^) 


8±1 


Angular radius of the shell R^ol (') 


17.5±3.0 


Linear radius of the shell (pc) 


18±5 


H2 mean column density (cm~2) 


(2.9±1.3)xl020 


H2 mass of the shell (M©) 


(8.2±1.7)xl02 



et al. 2004), and assuming a previous O-type phase with 
values of M = 2.3x10"'' M© yr"^ and v^c = 1950 km s"^ 
(Smith et al. 2002), we obtain a mechanical luminosity for 
the stellar wind of WR 157 during the O and WR phases, 
Lo ~3.5xl0^^ erg/s and Lw ~1.3xl0^^ erg s"\ respec- 
tively. Bearing in mind a lifetime to ~ 5x10^' yr for the O- 
type phase of a star with an initial mass of 40 Mq , and twR 
= 0.4x10® yr for the WR phase (Meynet et al. 1994), the 
mechanical energy injected into the interstellar medium in 
each evolutionary phase is -Ewq — 5.5 x 10'*^ erg and E^,^^ ~ 
1.6x10^° erg. 

The ratio e = E\^IE^ estimated considering only the 
WR phase of WR 157 is 6x 10"^. We can conclude that the 
stellar wind of the WR star alone is the main one responsible 
for shaping the interstellar bubble around the open cluster. 

The distributions of the interstellar dust and the ion- 
ized and molecular material in the environs of Mrk 50 can 
be explained as the consequence of the action of the stel- 
lar winds of WR 157. As regards Hi gas, no neutral atomic 
counterpart of SG 13 could be identified from the present 
study, since the spatial distribution of Shells A and B is 
very different from that of SG 13. It is likely that the small 
amount of Hi gas resulting from the photodissociation of the 
molecular gas would be hardly detectable against the strong 
background Hi emission. 

As regards the ionized gas, we can suggest two possible 
scenarios: 
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• The inner bright optical filaments were probably gener- 
ated by the action of the present WR phase of the star. The 
shock fronts have ionized and dissociated the molecular cir- 
cumstellar environment creating photodissociated regions. 
The outer optical filaments may have been swept up by the 
stellar winds in earlier stellar phases. 

• The other possibility is that the distribution of the dif- 
ferent optical filaments is far from being coplanar. If this 
were the case, we would be observing the projection of the 
filaments over the plane of the sky, the real dimensions of 
the structures being then larger than indicated in Table [2l 



7 CONCLUSION 

Using the CGPS high-resolution radio continuum and 21 cm 
Hi line data, supplemented by previous optical, IR and CO 
surveys, we arrive at the following conclusions concerning 
the ring nebula SG13 and the associated star WR 157. 

The radio continuum emission correlates extremely well 
with the optical DSS R and IR 60 /im images. Interference 
effects from nearby Cas A limit our ability to reliably de- 
termine the spectral index between the two CGPS radio 
frequencies, however the spectral index between 1420 MHz 
(CGPS) and 2700 MHz (Fiirst et al. 1990) is consistent with 
thermal emission. 

A partial ring of CO emission in the range -56 to -43 
km s~^ is seen to circumscribe the optical, IR (60 fim) and 
radio continuum emission of SG13. This range of velocities 
is consistent with previously measured velocity determina- 
tions by Lozinskaya et al. (1986) and Pedlar (1980), based 
on optical lines and the H166 radio recombination line, re- 
spectively. 

An analysis of the spatial distribution of IR point 
sources having colours characteristic of YSOs shows that 
an excess of such sources appears projected over the molec- 
ular ring surrounding SG13, suggesting that star formation 
triggered by the stellar members of Mrk 50, in particular 
WR 157, is taking place. 

The kinematics and dynamics of the gas, dust and 
molecular material around SG13 are entirely compatible 
with the hypothesis that the WR star alone is responsible 
for shaping the ISM around the open cluster Mrk 50. 

A 50' diameter Hi shell (shell A) is detected in the - 
32 to -50 km s~^ range to the north of SG13. Although 
this range partially overlaps the velocity range of the CO 
ring surrounding SG13, this Hi structure is entirely distinct 
from the other structures (molecular ring, optical, IR and 
radio continuum). The star WR157 appears projected onto 
the southern boundary of this Hi shell. An expansion ve- 
locity of 13±2 km s~^ is inferred from a velocity-position 
diagram. A second, smaller. Hi cavity (shell B) is also de- 
tected in the velocity range -52 to -41 km s~^, suggesting 
a lower expansion velocity of 9±2 km s~^. It too does not 
coincide with any of the structures associated with SG13. 
We conclude that these two cavities or shells are unlikely to 
be physically related to SG13. 

The lack of an Hi shell, which has been detected associ- 
ated with a large member of interstellar bubbles, is probably 
due to a low column density of the Hi gas resulting from the 
photodissociation of the molecular gas. 
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